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By using a suitable set of the surface energy coefficient, nuclear radius, and universal function, the 
original proximity potential 1977 is modified. The overestimate of the data by 4% reported in the 
literature is significantly reduced. Our modified proximity potential reproduces the experimental 
data nicely compared to its older versions. 
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I. INTRODUCTION 

o ■ 

<n : 

Recently, great theoretical and experimental efforts are taken to studying the fusion of heavy nuclei leading to 
several new phenomena including the understanding of the formation of neutron -rich and super heavy elements [J 0] ■ 
The precise knowledge of the interaction potential between two nuclei is a difficult task and continuing efforts are 
needed in this direction. This problem has been of very active research over the last three decades and remains one 
of the most widely studied subject in low-energy heavy- ion physics [H-Hoj. 

i 1 The total interaction potential is sum of the long range Coulomb repulsive force and short range nuclear attractive 

force. The Coulomb part of the interaction potential is well-known, whereas nuclear part is not clearly understood. A 
large number of efforts have been made to giving simple and accurate forms of the nuclear interaction potentials [J- 
• Among such efforts, proximity potential is well known for its simplicity and numerous applications. Based upon 

the proximity force theorem [2, [f| , a simple formula for ion- ion interaction potential as a function of the separation 
£h . between the surfaces of two approaching nuclei was presented [H, H| ■ 

As pointed out by many authors [7[, original form of the proximity potential 1977 overestimates the experimental 
H ■ data by 4% for fusion barrier heights. In a recent study involving the comparison of 16 proximity potentials, one 
■ of us and collaborators pointed out that proximity potential 1977 overestimates the experimental data by 6.7% for 
| symmetric colliding nuclei |IJ . Similar results were obtained for asymmetric colliding nuclei [l[ . 

With the passage of time, several improvement/ modifications were made over the original proximity potential 1977 
to remove the gray part of the potential. It includes either the better form of the surface energy coefficient @ or 
the universal function and/or nuclear radius [7]. A careful look reveals that these modifications/improvements are 
not able to explain the experimental data [ETIll- A deep survey also pointed out that these technical parameters 
(i.e. surface energy coefficient, nuclear radius, and universal function) were chosen quite arbitrarily in the literature. 
Among them, the surface energy coefficient is available in a large variety of forms from time to time [3, [2] . It affects 
the fusion barrier heights and cross sections significantly [H, 0. Also, nuclear radius is available in large variety of 
forms (J, |2| ■ These forms varies either in terms of its coefficients or either different mass or isospin dependence. The 
third technical parameter i.e, the universal function, is also parametrized in different forms [l], HUB], @] ■ Unfortunately, 
no systematic study is available in the literature, where one can explore the role of these technical parameters in 
fusion barrier positions, heights, and cross sections. Alternatively, a best set of the above-mentioned parameters is 
still missing. 

In the present study, our aim is to modify the original proximity potential 1977 by using a suitable set of the 
above-stated technical parameters available in the literature. In addition, to compare the final outcome with the 
huge amount of experimental data available since last three decades. The choice of the potential and its form to be 
adopted is one of the most challenging task when one wants to compare the experimental data with theory. The 
present systematic study includes the reactions with combine mass between A = 19 and A = 294 units. In total, 390 
experimentally studied reactions with symmetric as well as asymmetric colliding partners are taken into consideration. 
Section [TT] describes the Model in brief, Section [TTT] depicts the Results and Summary is presented in Section IIV1 
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II. THE MODEL 



The total ion- ion interaction potential Vr{r) between two colliding nuclei with charges Z\ and Z2, center separation 
r, and density distribution assumed spherical, and frozen, is approximated as Q 

V T (r) = V N (r) + ^^, (1) 
r 

where e is the charge unit. The above form of the Coulomb potential is suitable when two approaching nuclei are well 
separated. The nuclear part of the potential Vjv(r) is calculated in the framework of the proximity potential 1977 [1] 

as 

V N (r) = 47rfl 7 6$( - *~ 2 ) MeV, (2) 

where R = * s ^ e recmce< i radius. Here Cj denotes the matter radius and is calculated using relation Q 

Ci = Ci + -^U (i = 1,2), (3) 

where a denotes the half-density radii of the charge distribution and ti is the neutron skin of the nucleus. To calculate 
Ci, we used the relation given in Ref. [7J as 

/ 7 b 2 49 b 4 \ 
c i = R 00i l-^-Tnr + (* = 1,2). (4) 

Here, i?oo is the nuclear charge radius read as 

^ = 1.2332^ {l + ^ -0.151541 (^)) ^ (5) 

where Ni and Zi refer to neutron and proton contents of target /projectile nuclei. This form of radius is taken from 
the recent work of Royer and Rousseau [ll[ and is obtained by analyzing as many as 2027 masses with N, Z > 8 and 
a mass uncertainty < 150 keV. The neutron skin t{ used in Eq. © is calculated according to Ref. Q. 
The surface energy coefficient 7 was taken from the work of Myers and Swiatecki [l2j and has the form 



7 = 7o 



1-fc, 



N - Z s 1 
A 



(6) 



where N and Z refer to the total neutrons and protons content. It is clear from Eqs. ([5]) and ([6]) that both nuclear 
radius as well as surface energy coefficient depend on the relative neutron excess. In the above formula, 70 is 
the surface energy constant and k s is the surface-asymmetry constant. Both constants were first parameterized by 
Myers and Swiatecki [l2| by fitting the experimental binding energies. The first set of these constants yielded values 
70 and k s = 1.01734 MeV/fm 2 and 1.79, respectively. In original proximity version, 70 and k s were taken to be 
0.9517 MeV/fm 2 and 1.7826 [l3j], respectively. Later on, these values were revised in a large variety of forms depending 
upon the advancement in the theory as well in experiments [H 0- in total, 14 such coefficients are highlighted in 
Ref. Q and the role of extreme 4 sets is analyzed deeply. Out of them, two best sets of surface energy coefficients 
are stressed. In the present study, we shall restrict to the latest set of 7 values i.e. 70 =1.25284 MeV/fm 2 and k s = 
2.345 presented in Ref @- This particular set of values were obtained directly from a least-squares adjustment to the 
ground-state masses of 1654 nuclei ranging from 16 to 263 106 and fission-barrier heights [1J|. 

The universal function $( r ~ c j ) ~ c ' 2 ) used in Eq. ([T]) has been derived by several authors in different forms [3, [H, 0]- 
In original proximity potential, §{ r ~ Cl b ~ C2 ) was parametrized in the cubic-exponential form [|| 

f -3 (£ - 2 - 54 ) 2 - °- 0852 K - 2 - 54 ) 3 . fOT £ < L2511 . 

(e) ~ I -3.437 exp (-^) , for £ > 1.2511 , ( } 

with £ = (r — C\ — C-i)jb. The surface width b (i.e. b — -^a with a = 0.55 fm) has been evaluated close to unity. We 
labeled this universal function as $-1977. 
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Later on, Blocki et al., 5] modified the above form as 

{-1.7817 + 0.9270$ + 0.143£ 2 - 0.09£ 3 , for £ < 0.0 , 
-1.7817 + 0.9270$ + 0.01696£ 2 - 0.05148<f, 
for 0.0 < £ < 1-9475 , ( 8 ) 
-4.41 cxp (-g^Ts) , for £ > 1.9475 . 

In the present study, we use this form of universal function and marked it as $-1981. By using the above stated 
parameters, we construct a new proximity potential and labeled as Prox 2010. Along with the above modified form, 
we shall also use the original proximity potential 1977 [H and its recently modified form proximity potential 2000 0]. 
We labeled them as Prox 1977 and Prox 2000, respectively. 



III. RESULTS AND DISCUSSIONS 



By using the above new version of the proximity potential (Prox 2010) along with its older versions (i.e. Prox 1977 
and Prox 2000), fusion barriers are calculated for 390 reactions by using the conditions: 

dV T (r), n , d 2 V T (r), 

—^\ r =R B =0, and / 2 V ; \ r =R B < 0- (9) 
ar ar 

The height of the barrier and position is marked, respectively, as Vb and Rb- 

As one see from the preceding section, three factors govern the success of proximity potential are (i) the surface 
energy coefficient, (ii) the universal function, and (iii) nuclear radius. We analyzed the literature very carefully and 
found that the latest information on these three factors can shape the new proximity potential. Recently, the role 
of surface energy coefficient stated above is studied in detail in Ref Q. As for as radius is concern, we shall restrict 
to its latest form given in Ref. jllj . However, the role of the third parameter i.e., the universal function in fusion 
barriers is analyzed in Fig. 1. Here, we display AVb (%) and ARb (%) defined as 

Ytheor _ yexpt 

^V B (%) - ^—^^ x 100, (10) 



and 



ntheor _ r>expt 

Ai?B (%) = " B x 100, (11) 

as a function of Z1Z2 using two sets of above mentioned universal functions [Eqs. (0 and l[8])]. It is clear from the 
figure that deviations are significantly reduced by using $-1981 compared to its original form $-1977. The universal 
function $-1981 reduces the average deviation over 390 reactions by 1 % for fusion barriers. The experimental values 
are taken directly from the literature [H, [H, 0] . Actually, it is clear from the literature that no experiment can extract 
information about the fusion barriers directly. All experiments measure the fusion differential cross sections and then 
with the help of a theoretical model, one can extract the fusion barriers. 

In Fig. 2, we display the theoretical fusion barrier heights V| heor (MeV) and positions R t B eor (fm) verses the 
corresponding experimental values. We note from the figure that Prox 2010 potential reproduces the experimental 
fusion barrier heights within 1.4%. This result is in close agreement with other recently parametrized potentials 
presented in Ref. [l|. However, the original form of the proximity potential 1977 presented in Ref. [l| overestimates 
the data by 6.7% for symmetric colliding nuclei. However, the fusion barrier positions show some scattering from the 
central line (marked by shaded area). This scattering may be due to the variation in the experimental setups and 
theoretical method one used to extract these values 



ng may 

mm- 

We quantify our outcome in Figs. 3 and 4. In Fig. 3, the percentage deviations between the theoretical and 
experimental values are presented. The original proximity potential 1977 (Prox 1977) along with its recently modified 
form (Prox 2000) are also displayed. We note from the upper panel of Fig. 3 that Prox 2010 potential on average 
gives better results compared to its older versions for fusion barrier heights. However, slight deviations are visible for 
fusion barrier positions. This may be due to the fact that in the proximity potential Prox 2010, we use the value of 
surface energy coefficient that gives stronger attraction compared to one used in Prox 1977 and Prox 2000 potentials. 
Therefore, in Prox 2010 potential, the counterbalance between the repulsive Coulomb and attractive nuclear part of the 
interaction potential occurs at larger distances, and hence it pushes the barrier outwards. The fusion barrier heights 
are reproduced within ± 5% on average. On the other hand, fusion barrier positions reproduced the experimental 
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values within ± 10%. Especially for the heavier colliding nuclei, we see that Prox 2010 potential reproduces the data 
much better on the average compared to other versions. For lighter nuclei, however, small scattering is visible. This 
could also be due to the uncertainty in the radius of the lighter colliding nuclei. In Figs. 1-4, only 155 reactions are 
displayed to maintain the clarity. The average deviation for the fusion barrier heights over 390 reactions is 0.77 % 
using our modified potential Prox 2010, whereas Prox 1977, and Prox 2000 give 3.99 %, and 4.45 %, respectively. 
This shows that our modified proximity explains the experimental data nicely. 

In Fig. 4, we display the difference between the theoretical and experimentally extracted fusion barriers. We further 
note that Prox 2010 potential gives better results. The difference especially for the heavy systems is significantly 
reduced. This was the problem with original as well as its recently modified form as pointed out by several authors 0, 
[§J]. It is clear from Figs. 3 and 4, that Prox 2010 potential is able to reproduce the experimental data much better 
than its older versions. The small difference is not significant because of the uncertainties in the analysis of the 
experimental data. 

Finally, we test our newly modified proximity potential Prox 2010 on fusion probabilities. In Fig. 5, we display 
the fusion cross sections <Tf us (in mb) as a function of the center-of-mass energy E cm (MeV) for the reactions of 
2e Mg + 30 Si [l7|], le O + 46 Ti O, 48 Ca + 48 Ca [H, 12 C + 92 Zr M^Ca + 58 Ni $\, and le O + 144 Sm H- The 
fusion cross sections are calculated using well known Wong model [23[. The older versions of proximity potentials 
that is, Prox 1977 and Prox 2000 are also displayed. It is clearly visible from the figure that Prox 2010 potential 
is in good agreement, whereas, its older forms are far from the experimental data. We further note that Prox 1977 
and Prox 2000 potentials show similar results. It means that no improvements is seen in Prox 2000 potential as was 
claimed in Ref. [7]. 



IV. SUMMARY 



In the present study, we present a best set of the surface energy coefficient, the nuclear radius, and the universal 
function available in the literature. We find that these parameters which were used quite arbitrarily in past years 
affect the fusion barrier heights, positions, and cross sections significantly. By using the above set of parameters, a 
new proximity potential is constructed. Our newly constructed proximity potential Prox 2010 reproduces the fusion 
barriers and cross sections better than its earlier versions. 
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FIG. 1: The percentage deviation AVb (%) and ARb (%) as a function of Z1Z2 using two different sets of universal functions 
[Eqs. ID and ([5)1] implemented in the original proximity potential Prox 1977. The experimental values are taken from 
Refs. 
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FIG. 2: The fusion barrier heights Vb (MeV) and positions barriers Rb (fm) as a function of the corresponding experimental 
values using our modified proximity potentials Prox 2010. The experimental values are taken from Refs. [HHP 
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FIG. 3: The same as Fig.l, but for different older proximity potentials along with our modified form i.e., Prox 1977, Prox 2000, 
and Prox 2010, respectively. 
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FIG. 4: The variation of AVs (= V B heor 

2000, and Prox 2010 potentials. The experimental values are taken from Refs. [J HQ. 
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FIG. 5: (Color online) The fusion cross sections af ua (mb) as a function of center-of-mass energy E c . m . using older versions 
of proximity potential (Prox 1977 and Prox 2000) along with new version (Prox 2010). The experimental data are taken from 
Morsad 1990 \rft, Neto 1990 Stefanini \v^, Newton 2001 [U, Sikora 1979 [ij, and Leigh 1995 



